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This paper presents an analytical study of synchronization in an array of output- coupled temporal Boolean 
networks. A temporal Boolean network (TBN) is a logical dynamic system developed to model Boolean 
networks with regulatory delays. Both state delay and output delay are considered, and these two delays are 
assumed to be different. By referring to the algebraic representations of logical dynamics and using the 
semi-tensor product of matrices, the output-coupled TBNs are firstly converted into a discrete-time 
algebraic evolution system, and then the relationship between the states of coupled TBNs and the initial state 
sequence is obtained. Then, some necessary and sufficient conditions are derived for the synchronization of 
an array of TBNs with an arbitrary given initial state sequence. Two numerical examples including one 
epigenetic model are finally given to illustrate the obtained results. 

A central focus of modern biology concerns understanding the function of a cell, the manner how a cell 
performs its function and the ways in which cellular systems fail in disease. With the rapid development of 
systems biology' "', genetic regulatory networks have recently become a newly developing area and have 
attracted much attention due to its close relationship with neural networks, protein webs and other biological 
systems'"^. A genetic regulatory network contains a group of DNA elements in a cell which interact with others 
directly or indirectly through their mRNA or protein products. Thus, a genetic regulatory network governs the 
expression levels of the mRNA and the protein products. With time going, various types of genetic regulatory 
networks have been proposed, such as Markov-type genetic networks" and Boolean networks'. Each of these 
models has its own advantages and limitations in terms of the following considerations: the accuracy of the 
approach relate to the level of genetic, biological phenomenon being modelled; the experimental feasibility of 
model construction; and the computational complexity of the model inference from the available data. 

The concept of Boolean network (BN) was firstly proposed by Kauffman' as a simplified model of the behavior 
of large networks of randomly interconnected binary (on-off) genes. One of the main advantages of Boolean 
networks is the comprehensibUity owing to the simplicity of the representation. A BN can effectively model a 
genetic regulatory network with its gene products (the outputs) and the substances from outside or the envir- 
onment that can affect the gene (the inputs). In a BN, the state of each node has only two unique states: 1 (on) or 0 
(off). For a gene, "on" corresponds to the gene being expressed; For inputs and outputs, "on" corresponds to the 
relative substances being present. Moreover, each node evolves its state according to a Boolean function, which is 
a logical function, at each discrete time point. 

Since the BNs can provide a general description of the behavior of many living organisms, they have attracted 
great attention from researchers in different fields. There are various examples of genetic regulatory networks 
modeled using BNs, such as the yeast cell-cycle networks'; the control of the mammalian cell cycle'"; the yeast 
transcriptional networks"; the network predicting the expression of the segment polarity genes in the fly 
DrosophUa melanogaster'^, and so on. Recently, a new method called semi-tensor product (STP) of matrices 
has been firstly proposed by Cheng and his colleagues'''''''. By using the method of STP, a Boolean function can be 
uniquely expressed as an algebraic expression, and hence a BN can be expressed as a standard discrete-time 
algebraic system. This method has been proved to be quite useful to study BNs concerning the topological 
structure, such as the fixed points, cycles, basin of attractors and the transient time. Based on this approach, 
many interesting results have been obtained in the last few years, such as stabilization of BNs with impulsive 
effects''\ optimal control of logical control networks", etc. 

A curious ability of some real- world systems is that they can evolve in perfect synchronization. In the past 
decades, the phenomenon of collective behavior has inspired a large amount of research, such as synchronization 
analysis and control of complex networks'^ consensus in multi-agent systems^^'^^ synchronization of 
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Kauffman networks^'*, and cooperation of networks''^"^', etc. It was 
demonstrated that the nodes in a second-order Kuramoto model 
perform a cascade of transitions toward a synchronous macroscopic 
state"'*. A time delay, which is a source of instability, exists in various 
engineering systems such as chemical processes, or long transmis- 
sion lines in pneumatic systems. Some interesting results have also 
been obtained concerning the effect of time-delay on network syn- 
chronization'^'". A processing coupling- delay is effectively used to 
annihilate the quenching of oscillation in a network of coupled oscil- 
lators by switching the stability of amplitude death'-^. Recently, inter- 
ests have also been extended to synchronization of BNs, due to 
their potential applications in biology, chemistry, as well as 
engineering^''^*. 

Synchronization and controllability are structural properties of a 
system, and they are fundamental concepts in control theory. 
Recently, synchronization of two delay-free deterministic BNs and 
synchronization in an array of coupled BNs have been respectively 
studied in" " and the synchronization of multi-valued logical net- 
works have been studied in'^. In"', we investigated the synchroniza- 
tion in an array of BNs with time delay. It is well known that time 
delay is very common in real world and also is inevitable in biological 
systems and engineering problems. Moreover, in many situations, 
time delay can not be ignored in practice, since it may heavily affect 
the dynamics of the networks""*. One interesting kind of BNs, called 
temporal Boolean networks (TBNs) were developed to model regu- 
latory delays in''. TBNs own a more complex structure when the time 
delays in Boolean variables are allowed to be different. TBNs can 
approximate real cellular regulatory networks or genetic networks 
better than BNs. Hence, it is desirable to study synchronization of 
TBNs. In the past few years, many results about TBNs have been 
obtained, such as output controllability of temporal Boolean control 
networks (TBCNs) and complete synchronization of TBNs*". 
However, there has been no results investigating synchronization 
in an array of output-coupled TBNs, to the best of our knowledge. 
In*", the authors studied the synchronization of two coupled TBNs in 
a drive-response configuration. Sometimes, in the real world, we 
need to judge whether an array of networks can be synchronized 
or not. The output coupling phenomenon is ubiquitous in many 
biological networks and physical systems. For instance, a cellular 
system contains many different components to synthesis macromo- 
lecular in cells. In an idealized model of a metabolic feed-back control 
cycle, each unit component consists of three parts: a genetic locus, a 
ribosome and a cellular locus. Each of these three parts interacts with 
each other by its corresponding products, such as mRNA, protein 
and metabolic species. Furthermore, different components interacts 
with each other by its corresponding metabolic products, which is a 
type of output coupling. Thus, investigating synchronization in an 
array of output-coupled TBNs is meaningful and challenging, and 
can be helpful to disclose the control mechanisms for macromole- 
cular synthesis and the behavior of a cell. 

Motivated by the above discussions, in this paper, we focus on 
synchronization in an array of output-coupled TBNs. Our main 
objective is to obtain a general analytical approach to study syn- 
chronization in output-coupled TBNs. Both state delay and output 
delay are considered in the output-coupled TBNs. The main results 
are derived by using the method of semi-tensor product of matrices 
which can convert a deterministic BN into a unique algebraic rep- 
resentation. By using the algebraic representations of TBNs, we 
derive some necessary and sufficient conditions for the synchroniza- 
tion in an array of TBNs. 

Some basic notations: The standard notations will be used in 
this paper. Throughout this paper, the notation Ij- denotes the k- 
dimensional column vector with all entries equalling to 1 . Denote the 
set of two logical variables by 27 = { 1 ,0} , where 1 ~ T means "true" 
and 0 ~ P means "false". A logical variable u will take a value 
from D, which is expressed as u eD. Define a delta set as 



A„ : = {(5^j! = 1,2, ...,«}, where is the fth column of the identity 
matrix I„. We denote the fth column of matrix A by Col,(A) and 
denote the set of columns of matrix A by Col(A). We also denote the 
fth row of matrix A by Row,(A) and denote the set of rows of matrix A 
by Row(A). An n X m matrix L is called a logical matrix if the 
columns of L, denoted by Col(L), are in the form of 5*. That is to 
say, Col(L) ^A„. Denote the setofn X s logical matrices by£„xs-If 
LeL„ X s, it can be expressed as L = [S'^ ,5'^ , . . . ,5'^]. For simplicity, it 
can be denoted as L = S„[ii, z'J. We use ®^[X,(t) to denote 
Xi(f)®X2(0® ■ ■ ■ ®Xm(0. and use K^i7,(0 to denote 
yi{t)Ky2{t)K-- -KyNit). 

Results and Methods 

In this paper, we study the following array of M output-coupled BNs, 
with each one being an N-node system: 

' x'.{t+i) =f;{xl{t), . . . ,xf (0,x,i(t- 1), . . . ,xf{t-i), 

Xl{t-z,),...,Xf{t-n),y,{t),---Mt)), 
i=l,2,...,N 

y,it) ^gj[xj{t), . . . ,Xf{t),X]{t- 1), . . . ,Xf (f- 1), . . . , 

X]{t-X2 + l),...,Xf{t-X2 + l),X]{t-Z2),..., 

Xf(t-T2)),; = 1,2,...,M 

where is the fth node of the jth TBN, y/ is the binary output, 
// : {l,0}'^' + "^+^^{l,0}, gj : {1,0}<^^ + ''''^{1,0} are Boolean 
functions, and Ti, T2 are positive integers describing the state delay 
and the output delay. In many biological networks and genetic 
regulatory networks, time delay is ubiquitous due to finite signal 
propagation time, finite chemical reaction times, memory effects, 
and so on. The time delay Ti describes delay occurring between states, 
while time delay T2 describes delay resulting from output. For 
example, in a genetic regulatory networks, the level of a gene maybe 
affected by its past n gene levels. Different genes interact by its cor- 
responding gene products, which may cost time to be transmitted 
from one gene to another one. Thus, model (1) can be an efficient 
model to describe some biological networks or genetic regulatory 
networks with state delay and output delay. We simply denote 

Xj{t) = (xl{t),Xf{t),...,Xf{t)y the states of the TBNs ( 1 ) at time 

instant t. We observe that the state evolution of the output-coupled 
TBN (1) depends on the initial state sequences, and there are two 
cases as follows. Case (I): s Z2, the initial state sequence is X^(0), 
X/-1), Xj(-Ti) for; = 1, 2, M; case (II): Ti < T2. the initial 
state sequence is Xj{0), X,(- 1), . . ., X,(-T2) for; = 1, 2, . . ., M. Now, 
we define complete synchronization in the array of TBNs (1) as 
follows: 

Definition 1 The array of TBNs are said to be synchronized if for 
every initial state sequence, there exists a positive integer k, such that t 
> k implies Xjit) = Xj{t) for every distinct 1 ^ i,j ^ M. 

Remark 1 In this definition, the value of k depends on the initial 

state sequences. Nevertheless, since the {1,0}^^' ^^''^ (or {1,0}'^^* ''^J 
is a finite set, we can always choose a constant k large enough which is 
independent of the initial states sequence. 

In the following, we establish some necessary and sufficient cri- 
teria for synchronization in an array of output-coupled TBNs. To 
acquire a strict analysis, we ought to apply algebraic representations 
of the TBNs according to the STP theory. 

Let Xj{t)=X]{t)KXj{ty---KXf{t) for ; = 1, M and 
y[t) = yi{t)Ky2{t)K ■ ■ ■KjKAf(t). Let F' and Gj be the structure matrix 
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of fl and gj. Then we can make Aj{t,T2) = 
Xj{t)Xj{t-l) ■ - -Xjit-Ti) forj = 1,2, ...,M to obtain the algebraic 
representations of y{t). By equation (1) and the algebraic representa- 
tions of Boolean functions, we have 

yj{t) = GjXj{t)Xj{t-l) ■ ■ ■Xj{t-Z2) = GjAj{t,Z2). (2) 

Hence, we get y{t) = (^®f^^Gj>j>f^,Aj{t,Z2) = GKf^iAj{t,T2), 
where G= To obtain the accurate algebraic expressions of 

><jliAj(t,T2), we also make M;(tT2)=Xj(f-l) • • •X/f-T2) to facil- 
itate the expression. Then we have 

Kf^^Aj{t)=X,{t)Xi{t-l) ■ ■ ■Xi(f-T2)X2(0 ■ ■ ■X2(t-T2) 

■■■XM{t)---XM{t-T2) 
= Xi (OMi (t,T2)X2 (OM2 (t,T2)X3 (OM3 (t,T2) 

■■■XM{t)MM{t,Z2) 
= Xi(0W[2N,2'2«lX2(0Mi(t,T2)M2(t,T2)X3(0 

M3{t,Z2)---XM{t)MM{t,T2) 
= (/2« ® W|2N 2'2A'])Xi(0X2(0Mi (t,T2)M2(f,T2) 

X3(0M3(f,T2) ■ ■ ■XM(0MM(t,T2) (3) 

= xflT ' (^2"V ® W[2A.,2''2«| ) (x*! iX^(0) Ml (t,T2) 

M2{t,Z2)---MM{t,Z2) 
= Wo(K^lX^(0)Mi(t,T2)M2(t,T2)---MM(t,T2) 

= Wo(KjliX;(0)Xi(t-l) ■ ■ ■Xi(t-T2) ■ ■ ■XM(t-l) 
•■■XM(f-T2), 

where Wq = xfl^ ^ (I2'N ® 2i'2«] ) • By induction, we yield 
Kf^,Aj{t,T2) ^ mKfL^Xj{t)Kf^,Xj{t-l) . . .KfL^X,{t-Z2), (4) 

where Wi= WoXjLT Xh-'^-" ® W,) , = '^/2'N ® W^^, ^,,(.^_,)^^ , j 

= 1, T2 - 1. Hence, we have y{t) ^ GWiKf^^Xj{t)>ji^ 
Xj{t-l)---xfL,X,{t-Z2). 

Now, we try to obtain the algebraic representations of the states 
Xj(t + 1). By equation (1) and algebraic representations of Boolean 
functions, we have Xj(t-|- 1) =FjX;(t)Xj(t- 1) ■ ■ ■ Xj{t-zi)y{t). By 
denoting C(t,Ti)=Xj(t)Xy(t-l) ■ ■ ■ Xj{t-Zi)y{t) for simplicity, we 
then get 

Xj{t+l)=FlC{t,Zi)FjC{t,Zi) ■ ■ ■ FfC{t,z{) 

= f;^(/2(,,.,),v+«®f;')<l'(x, + i)N+MC(t,Ti)---PfC(Mi)(5) 
=Ff[Uz^)=F,Xj{t)Xj{t-l) ■ ..X,{t-z,)y{t). 

where Fj = -F/Kjl2 { {^2^^^ + '^(^1 + i)n+m | ■ 

That is to say, we have obtained the following equivalent algebraic 
representations of output-coupled TBNs (1): 

x,(( + 1) =F,x,(Ox,(t- 1) . . .x,(t-TiMO 

y{t) = GWiKf^ ,X,[tyf^ ,X,(t - 1) ■ ■ -Kf. 1 X/t - Z2) ^ ^ 



where Fj is an a'^ x 2(^' + i)'^+*' matrk, G is an 2** x 2(^^ + 1)'^** 
matrbc and Wi is an 2(^^ + 1)^** x a'^^ + D'^** matrix. 

From equation (6), it is easy to see that the state of each TBN is 
determined by its delayed states and the output y{t), which means 
these TBNs are coupled through the output y{t). Let 
KjliXi(t) =Xi(t)lxX2(0x ■ ■ ■kXm(0. using equation (6) and math- 
ematical iteration, we can derive the relationship between lxj*liXj(t) 
and corresponding initial state sequence. 

Before we derive the necessary and sufficient conditions for syn- 
chronization, we need the following proposition. From this proposi- 
tion, we can acquire the relationship between the state Xj{t) at time t 
and the initial state sequence which will be employed to obtain the 
criteria for synchronization of TBNs. 

Proposition 1 Let F^^f^fy W= W^^„_^(,,,,)«^ xfl^ 

|(^72'«®W[2M,2'('.+0"])'^m} ""'^ W2 = W'oK]l-\l2jMMlg)W)), 

where =xfl7i (^/2- ® W^^,, ^,(.,_,),^ ^ , 7 = 0, 1, ii - 1. 
Then we have 

KfLiXj{t+l)= FWGWi {l2{,, . ® W2) xjl iXj{t)x.fL iXj{t-l) 

■ ■ .><f^,XXt-T2)K^iX,(0xf.i^;(f-l) (7) 

■•■xf.lXXf-Tl). 

Remark 2 According to equation (7), we need to know whether Zj is 
larger than Z2 or not by considering these two cases: Case (I): Zi a Zj, 
which means that the initial state sequence isXj(O), Xj(—Zj)forj = 
1, 2, M; Case (II): Zj < Z2, which implies the initial state sequence is 
Xj(0), Xj(—Z2)forj = 1,2,..., M. In the following sequel, according 
to different case, we derive some necessary and sufficient conditions for 
synchronization in an array of output-coupled TBNs. We firstly con- 
sider Case (I), and then Case (II). 

Case I: Tj a T2. Now let us first consider Case (I): Zi s Zj, with the 
initial state Xj(0), Xj{ — Zi) for j = 1, 2, M. Since Ti and Z2 
describe the state delay and the output delay, in this case, it means the 
state delay is not smaller than the output delay. In a cellular system, it 
means the delay caused during the processes of the transcription of a 
gene is not smaller than that caused during the processes of the 
translation of a gene. Based on Proposition 1, we obtain the 
following proposition, which is derived to show the explicit 
relationship between the state Xj(t) and the initial states Xj(Q), 
Xj(-zi). 

Proposition 2 Lef ©1 =FWGWi (i2(,, + i)nm<E)W2^<^(,, + i)nm and 
El = 01 + <I>TiWM. then we have 

KfL iXj{t+l)= ©iKf^ ,X,(0k*1 iX,(f- 1) ■ ■ -xfL rXj{t-Zi) (8) 
and 

KfL ,m = ih-"-' ® ,mK-i^j{ - 1) • • -Kf. , X,( - Ti) }. (9) 

In the first case, i.e. Ti s which means the state delay is larger than 
the output delay, we can observe the value of states at time t -\- 1 
depends on its past Ti + 1 states. Then, using mathematical analysis 
and STP, we can obtain that the value of states at time t has a specific 
connection with the initial sequence, K^iXj{Q), K^iXj{ — l), . . . , 
iXj( — Ti). The left side of equation (9) is the STP of state of each 
TBN at the discrete time point t, while the right side is product of 
several matrices and the initial state sequences. Since (72mn ® lypVM) 
and Hi are constant matrices, we can present our necessary and 
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sufficient criterion for synchronization in the array of output- =K*liX,(t) ■ ■ -Kjli Xy(t- ti) to analyze the synchronization 

coupled TBNs as following. problem. According to equation (8), we obtain another form of rela- 

Theorem 1 Let equation (6) be the algebraic representations of the tionship between the state Xj{t) and the initial state sequence X,(0), 

array of TBNs (1). Then synchronization occurs iff there exists a x,<-l), X/-Ti) by mathematical induction. 

positive integer k satisfying 1 si k ^ kg such that Proposition 3 Let ©i defined as the same as in Lemma 2, then we 



( . _-,\ 1 gei inui 

Co/[fev« (x)1J,..«)h;1 £ d^, -.^ = 1+ ^ ii , ' ,i = U2, . . . .2" (10) „ ,n « 



where 



get that 

: ■ - - ■ ]\ : (12) 

where 



ko = min | ' > 1 ■ {h"-'-' ® 1 JnNM ) E\ = (/jnm ® lj.,,v,« ) H{ j > i| . (11) 

Remark 3 In equation (11) in Theorem 1, we can use 
k'o = minl^i>l : S,[—'E{,j>i^ to simpUfy it but meanwhile the 

restriction is strengthened. By k'g — min^i> 1 : S,\ =E{,j>i^, there 
must exist an integer j>k'Q such that Hj°=H-'j. Then we have 



2~-l 



lfilK;^,_,{[l,(.,.,)-.®(6>'lif')]%,+l)NM} (13) 
^|^2('i+^-0"*' 2('i + ')"■«] '^{^i+2-t)N M, 3 < t< Ti + 1, 

^ Lfli^r]:l^ { [/^f,, .,).v .« ® (©'lif ' ) ] 0(„ t > Ti + 2. 

(J2nm(x)1J,,„„)Hj° = (J2nm(x)iJ,,„h)E-'j. Hence, we can conclude By Proposition 3, we then present anothernecessary and sufficient syn- 

that ko < k'o- Equation (10) implies that after finite time steps, i.e. kg chronization criterion for the array of TBNs. 

steps, we can judge whether the array of TBNs can be synchronized or Theorem 2 The array of TBNs can be synchronized iff there exists a 

not. If this array of TBNs can not be synchronized before the kg-th positive integer k such that 
steps, then it can not be synchronized at any time. Thus, Theorem 1 

provides an efficient criterion for synchronization in an array of out- Qgi ( 0*i®i ^ c / ^\ ^ ■Xi = i-\- iinAli^ 
put-coupled TBNs. \ i k j ^ i"--' 

RemsLik 4 If the coupled TBNs can be synchronized after certain @ 

, , , , , , ^ ; , ., ^ where L, ts denned m Proposition 3. 

timesk,thenwehaveX\(t)=X'Jt)forl<=ij^i'<M,\<=ij^i'<=N „ ,. . -r , . .., 

J J ^ '■' •' Remarks Accoramg to equation (12), if we have an initial sequence, 

and t ^ k. Since the set of state values is a finite set, i.e. . m v /n\ m v i .\ M \ d 

r , , ,„i '-e- Kr=A(%f=i^;(-l)---I><r=i^i(-^i)=^2""('.")' 

< 521^^2"' ■ ■ ■ '^2" f ' suppose each synchronous state to be \ k @^ 

l. J obtain that the value of states at time k is Coli 0jL^"' . Therefore, 

Xj = d'2N, 1 < j < M, 1 < ! < N. Because of the definition o/ . _ -, I- J 

^^,Xj(t) =Xi (Ok ■ ■ ■KXM(t)> we can obtain the synchronous set of i®^^^ J " ^''^ total possible state values at time t In order to judge 

r . (j—i'j(2'^'^—i) 1 whether the coupled TBNs can be synchronized or not, we just need to 

Xj^iXj[t), which is <(52;Mw:Ai = H 2^ — 1 '' = 1'2, ...,2 >. know whether there exists a positive integer k such that the possible 

. J. ^ ^. /„! .r 1 • -x- I • state values are in the synchronous set 

According to equation (9), if we nave an initial sequence, i.e. i ■ \f mn \ ^ 

><^iX,.(0)><fiiX,(-l)---><f^iX,(-Ti) = 5;,,«(,„), we can obtain i^.^ : a, = 1 + ^ ~ = 1,2, . . . ,2^ I. 

KjliXj(0 = Co/,[(72,vA,(g)lJ„«„)H;]. ™«s, Co/[(/2N*.®lJ„«„)H'i] ^ 

is the total possible state value ofK^iXj{t). In order to judge whether Case II: Tj < Tj. Now we turn to case (II): Tj < T2. where the initial 

the coupled TBNs can be synchronized or not, we just need to calculate state sequence is X/0),X/-1), . ...X^-Tj) for; = 1,2, ...,M. In this 

whether there exists a positive integer k such that all possible state value case, it means the state delay is smaller than the output delay. The 

of K^iXy(A:) are in the synchronous set following proposition is firstly presented to explicitly express the 

j^^. _ jj^2A' N _ j^-j ^ relationship between the current state and the initial state sequence. 

5'^. -1-1+ ^^3^ ,i=l,2,...,2^|. Proposition4Let02=PWGWi(72(.,,ON«®M^^W'^2(n+O««>+O"^^ 

As we know. Hi is a logical matrix and so is the matrix of ^ ^^^^ ^ ^nd H2 =02Wr . then we have 

/2AIW ® ijti NM. Moreover, if (iii),-, = 1, then we have L ^ ' J 

(72««®ly,A'M)Col;(Hi)=Col,(/2''«®l2''.,N«) which means taking ^^^Xj{t + \) = &2<x!^^^Xj{t)^^^Xj{t-\) ■ ■ -^^^Xjit-Xi) (15) 
the !th column of matrix I2H .« ® 1 J, ,v m . Then, we can easily derive the 

following simplified criterion for synchronization in an array of out- ^^'^ 

put-coupled TBNs. <.,X,(0 = (4-.® i,';„M)Ha><)l,X,(0)><f^,X,(-l) . . .>^"^,X,(-t,)|. (16) 

!=1,2,...,2'^ ,;=!,..., 2"'^*'}. Then the array of output-coupled = 1.2, M. Thus, the STP of state of each TBN is determined by the 

TBNs (1) can be synchronized iff there exists a positive integer k such initial state sequence and a transition matrix (/2aj« ® ijtj ,vm)H2. Now 

that RoWi (Hj) =(0 0 0) fej 1 2 2^^' + i)n a( | /q^ we can present the following synchronization criterion which is similar 

^ V ' ' ' with Theorem 1. 

2(ti + i)NM Theorem 3 Let equation (6) be the algebraic representations of the 

In Proposition 2, we use the variable substitution array of TBNs (1). Then synchronization occurs iff there exists a 

Zi{t)=tK^^^Xj{t) ■ ■ -xfLj^Xjit-Ti) to analyze the problem of syn- positive integer k satisfying 1 < A: < fco such that 
chronization in an array of output-coupled TBNs. Based on the ^ (;_i)(2*'"-l) ) 

variable substitution, we derive the necessary and sufficient criterion Co;[(/2.v,«®i2.2»")^2] -'j ''J"" M = '^-\ — \ -'^i'^ — '2' r (1^) 

for synchronization, i.e. Theorem 1. In the following sequel, we use 

the mathematical induction without using the variable substitution where 



SCIENTIFIC REPORTS | 4:6292 | DOI: 1 0. 1 038/srep06292 



4 



Figure 1 | Boolean model of a coupled macromolecular synthesis in a cell with each unit consisting of a genetic locus (Ii or Lj), a cellular structure (Rj or 
R2) and a cellular locus (Ci or C2). Arrows correspond to positive interactions, circles to negative ones. 



to = mm I i : i > 1 , (/2.V *> ® ij,, .v «. ) Ej, = (/2» « 0 1^ » ,« ) E'j /or some j > i | . (18) 
Similar to Theorem 1, we can also obtain the same synchronous set, 



which is 



(f-l)(2* 



According to equation (16), we can obtain the total possible state 
values as Col [ (l2« m (8) 1 n m ) SJ^] ■ As long as there exists a positive 
integer k such that the possible state values at time k are in the 
synchronous set, then we can draw a conclusion that the coupled 
TBNs can be synchronized. 

Remark 6 In equation (18) in Theorem 3, we can also use 

k!a = min^i>\ :H2 = H^J>!| to simplify it but meanwhile the 

restriction is strengthened. And we can also get fco<fc'o. Equation 
(17) implies that after finite time steps, i.e. kg steps, we can judge 
whether the array of TBNs can be synchronized or not. If after fcg 
steps, this array of TBNs still can not be synchronized, then it can 
not be synchronized at any time. Thus, Theorem 3 also provides an 
efficient criterion for synchronization in an array of output-coupled 
TBNs. 

Since (-r2''*'® iJ'i'Vm) and 22 are logical matrices, we can derive 
the following simplified criterion for synchronization in an array of 
output-coupled TBNs. 



r f!-i)f2^*'-i) 

CoroUary2ietn2=|(A,— 1)2'^'^**+; ; ki = l+- ^^p^^ , 

i= 1,2, . . . ,2^ ,j=l, . . . ,2^' ^ Then the array of temporal Boolean 
networks can be synchronized iff there exists a positive integer k such 

thatRow,{^\) = (0,0, . . . ,0) , !e|l,2, . . . ,2(^^ + 1' ^ '^^ j Q.^. 

2(t2+l)NM 

Similar to Proposition 3, we can also use mathematical induction 
to obtain following proposition which implies the relationship 
between the current state and the initial state sequence by using 
Proposition 4. 

Proposition 5 Let ©2 be defined the same as in Proposition 4, then 
we can obtain that 

Kf^,Xj(t)=@\LfMf^,x,{oyf^,Xj{-i) . . .KfL,x,{-T2) (19) 

where 



Table 1 | Boolean functions of the Boolean 


model of coupled macro- 


molecular synthesis in cells 




variable 


Boolean function 


xl =genetic locus Li 


-^xlAy2 


jcj = ribosome Ri 




jCj = cellular locus Ci 




xl = genetic locus Li 


^xl Ayi 


= ribosome R2 


xl 


xl = cellular locus C2 





t = 2. 



if-lXj=(-2{ [l,(,,.>).M(S)(&2Lf')]<^U, + l)NM} (20) 
W^^(,^^,_,)HM 2{,2 + i)NM^^(i2+2-t)N 3<t<T2-|-l, 

, Lfiy,Ztli'{[l,(.^.).M^(&2Lf')] + l)NM},t>T2+2. 

According to Proposition 5, we can obtain following synchronization 
criterion to judge the array of output-coupled TBNs can be synchro- 
nized or not, which is similar to Theorem 2. 

Theorem 4 The array of TBNs can be synchronized iff there exists a 
positive integer k such that 

(!-l)(2"~- 



Co/ &ivi 



2N-1 



,! = 1,2, 



(21) 



where L®' is defined in Proposition 5. 

According to equation (19), we obtain the total possible state 

values are Col^Q^^-f"^- Following similar discussions in Remark 
5, the coupled TBNs can be synchronized if and only if there exists 
an integer k such that the possible state values are in the synchron- 
ization set. 



Numerical Simulation 

In this section, we present two numerical examples to illustrate our 
main theoretical results, which include a practical epigenetic 
example. 

Example 1 Let us first consider the BN which is used to model a 
molecular biological system*\ It is an epigenetic model proposed by 
Goodwin*'' of control mechanisms for the macromolecular synthesis in 
cells including: genetic locus Li, cellular structure Ri (ribosome) and 
cellular locus Ci, i = 1, 2. As is shown in Fig. 1, Li synthesises mRNA 
(Xj) to the cellular structure R,-, where its activity results in the syn- 
thesis of a protein, denoted by 7;. The protein then travels to some 
cellular locus Ci. A fraction of a metabolic species Mi generated by the 
protein returns to the genetic locus Li. Moreover, a reciprocal inter- 
action occurs from L2 to Li. Based on the analysis in", the variables 
and Boolean functions of the model are given in Table I. 

i'^f/l=(f////i>KA>'2,x;,xf),/2=(f2\/|,/|) = KA^,,X^,X^), 

^1 = —'X[, g2 = ^x'2- During the process of synthesis in cells, several envir- 
onmental conditions such as temperature, growth rate or concentration of 
nutrition, may cause a time delay of state and output Thus, it is reasonable 
to use TBNs inform (1) to model the network. Meanwhile, we can use the 
following case to illustrate the impacts of time delay occurred in the process 
of synthesis in cells towards the synchronization problem. Put 
xj{t)=xj{t)t><xj{t)t><x^{t) forj = 1, 2, and y{t)^yi{t)Ky2{t). In the 
TBNs, we let the initial state sequences be Xi(—l) = (1, 1, 1), X2(— 1) = (0, 
0, 0), xM = (0, 0, 0), X2(0) = (1, 1, 1). 

Consider: Xi = 12 = 1 which means that the state time delay is equal 
to the output time delay. 
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Figure 2 | (a): Time evolution of TBNs (22) with the given initial state sequences. The black sohd Une presents the state evolution ofxi, while the red dash 
line presents the state evolution of x^. (b): Total synchronization error £(t) of the coupled BNs (22) with the given initial state sequences. 

By using the STP of matrices, we also obtain the following algebraic nized by Theorem 1. Define the total synchronization error: 



(22) 



representation: 

{ xi{t+l)=Fxi{t-l)y2{t) 
X2{t+l)=Fx2{t-l)yi{t) 
y,{t) = Gx,{t-l) 
yy2(t) = Gx2{t-l) 

where F = 8f,[5, 5, 1, 5, 6, 6, 2, 6, 7, 7, 3, 7, 8, 8, 4, 8] and G = Ss[2, 1,2, 
1,2, 1,2, 1]. 

In order to see whether these two TBNs can be synchronized or not, 
we should calculate according to Theorem 1. Direct computation 
gives that Ej = H[^, then we have k'^ — S. According to Remark 3, we 
have ko<k'o = 5. Moreover, we obtain that (72« ® iJO'^i — 
{<5^4|! = 1,10,19,28,37,46,55,64}. Hence, the TBNs can be synchro- 



E{t) = l^'i (0 "^2(0 1 ■ P^8- ^('^) shows the time evolution of 

the TBNs and Fig. 2(b) shows the total synchronization error E{t). 
By Proposition 2, we also get Hi. Due to the complexity of this matrix, 
we just present the indexes of the first 50 columns as follows: 



Hi = .5212 [2305, 2305, 2369, 2369, 2433, 2433, 2497, 2497, 
65, 2433, 129, 2497, 193, 2817, 2817, 2881, 2881, 
2945, 2945, 3009, 3009, 2817, 513, 2881, 577, 2945 
641, 3009, 705, 3329, 3329, 3393, 3393, 3457, 3457, 
3521, 3521, 3329, 1025, 3393, 1089, 3457, 1153, 
3521, 1217, 3841, 3841, . . .] 



(23) 




(a) 



1000 2000 3000 
column index 



4000 





68 r 




64 




60 




56 




52 




48 




44 


X 


40 


0 




-a 


36 




32 


0 






28 ■ 




24 




20 




16 




12 




8 
4 







(b) 



1000 2000 3000 
column index 



4000 



Figure 3 | (a): The whole row indexes of each column of matrix Hi. (b): The whole row indexes of each column of matrix (726® lj6)Hi. Each point 
corresponds to the row index of each column, which implies the position of 1. 
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Figure 4 | Two coupled BNs with each network consisting of two nodes and one output. Arrows correspond to interactions between nodes and outputs. 



and Fig. 3(a) shows the whole row indexes of each column of matrix 
Hi, Fig. 3(b) shows the row indexes of each column of matrix 

(/26®l5)H=. 

Example 2 In Fig. 4, we have a system consisting of two coupled 
TBNs with delay Zi = = I, each TBN with two nodes. And its 
dynamics is described as follows: 

x\{t+l)=[x\{t)^x\{t)\ ^y2{t), 
xlit+l)=[x\{t)^xl{t)], 
yi{t)=x\{t)yxj{t-l). 
l{t+l)^x\{t-l)^yi{t). 



(24) 



xl{t+l)=xl{t)Axl{t-l), 
J2{t)=xj{t-l). 

Denote Xj{t) = x^^ {t)K>^ [t) for j = 1,2 and y{t) =yi{t)Ky2{t), then we 
can express the system into its algebraic form by using the STP of 
matrices as follows: 



and 



- x\(t+l)=F\xi{t)xi{t-l)y{t) 
xl{t+l)=F\x,{t)x,{t-l)y{t) 
yi{t) = GiXi(t)xi{t~l) 



- x\[t+l)=F\x2[t)x2[t-l)y[t) 
xl{t+l)=Flx2{t)x2{t-l)y{t) 

y2{t)=G2X2{t)x2{t-l) 



(25) 



(26) 



We obtain the explicit expression of the matrices F\,F\,F\,F2, Gi and 
G2, by using STP. Here, we omit the procedure of computation for 
simplicity. 

As we can see this system belongs to case (I). By Proposition 2, we 
obtain the matrix and Fig. 5 shows the row indexes of each columns 
of matrix Ei. In order to judge whether these two coupled TBNs can be 
synchronized or not, we should calculate E[ according to Theorem 1. 
Direct computation shows that (72" ® ij*) '^1 ~ (^2-' 0^5) S}' by using 



c 
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Figure 5 | The whole row indexes of each column of matrix Si obtained in system (24). Each point corresponds to the row index of each column, which 
implies the position of 1. 
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Figure 6 | (a): The whole row index of each column of matrix (72* ® 1 ) 5i obtained in system (24). (b): The whole row index of each column of matrix 
(724®lj4)Hj. Each point corresponds to the row index of each column, which implies the position of 1. 



the equation (11), then we have ko = 8. However, we get thatS.^ = S,\*, 
while Hj^H^. By using k'o = min^i>l : Hj=H{,j>i|, we have 

k'o = 9 which verifies Remark 3. It is easy to verify that 
Col[{l2>®\l,)'B.\] ^{,5^4 : ! = 1,6,11,16} /or 1 < f < 8. In Fig 6 
to Fig. 9, we plot the whole row index of columns of each matrix 
(^/2'' ® iJO^'i' ' ~ ^' 2' 8- F^ofn these figures, we can also draw 
the conclusion that Col[{l2i®l2i)S.[] ^{^i ■ ! = 1,6,11, 16}/or 1 < 
t s 8. Hence, this output-coupled TBNs can not he synchronized 
according to Theorem 1. 

In the following, we can also define the total synchronization error 
of the output-coupled TBNs: E{t) = \x\{t) -xl{t)\-\-\xl{t) - xl{t)\. 
Fig. 10(a) shows the time evolution of the state Xi(t) and X2{t) of the 
array of TBNs (24), and Fig. 10(b) shows the totaly synchronization 
error E(t) of the TBNs (24). From the state evolution and the total 
synchronization error of the TBNs (24), we can draw the conclusion 



that the TBNs can not synchronized. Our analytical results are well 
illustrated by this example. 

Remark 7 Due to the complexity of computation. Fig. 10(a) and 
Fig. 10(b) just plot time evolution and the total synchronization error 
of the TBNs with time up to t = 100. To some extent. Fig. 10(a) and 
Fig. 10(b) can illustrate that the coupled TBNs can not be synchro- 
nized. In fact, we have already made the computation of these TBNs 
with time up to t = 10000, observing that these TBNs can not be 
synchronized. Since the simulation result can not be well expressed 
in the figures with time t = 10000, we just plot Fig. 10(a) andFig. 10(b) 
till time f = 100 only. 

Conclusions 

In this paper, we have studied synchronization in an array of output- 
coupled TBNs. Both a state delay and an output delay Z2 are 
considered in the model. In a genetic regulatory network, state corre- 
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Figure 7 | (a): The whole row index of each column of matrix (724 ® 1 J4 ) Hj obtained in system (24). (b): The whole row index of each column of matrix 
(724®lj4)Hj. Each point corresponds to the row index of each column, which implies the position of 1. 
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Figure 8 | (a): The whole row index of each column of matrix (72* ® ijt ) obtained in system (24). (b): The whole row index of each column of matrix 
(72*®lj4)Hj. Each point corresponds to the row index of each column, which implies the position of 1. 



sponds to the level of gene, while the output corresponds to gene 
products. In real world, many evolutionary processes may experience 
external or internal changes. These changes maybe due to the envir- 
onmental disturbances or the interconnections of other genes, which 
can heavily affect the evolutional processes of genes, such as the 
transcription process or the translation process. By using the semi- 
tensor product of matrices, we first convert the TBNs into an unique 
discrete time algebraic system. Based on the equivalent algebraic 
representations of TBNs, we then study the coupled TBNs by con- 
sidering the two cases: Zi a T2 and Ti < T2. For each case, we derive 
necessary and sufficient synchronization criteria for the array of out- 
put-coupled TBNs. Different types of criteria are derived from dif- 
ferent aspects. These criteria provide efficient ways to judge whether 
an array of output-coupled TBNs can be synchronized or not. And 



they also implies that if the array of TBNs can not be synchronized 
among the upper bound time, then they can not be synchronized at 
any time, which means we only need to judge finite time steps. 
Finally, two illustrative examples including a practical epigenetic 
example are given to show the validity of our theoretical results. 
One interesting topic in the near future is to design logical controllers 
to make coupled TBNs synchronized. 

This mechanism is expected to be applicable to a wide range of 
biological systems, because synchronization phenomenon is very 
commonly in biological systems. Moreover, it is expected to help 
to understand the nature of cellular function, the behavior of genes 
and the manner in which cells execute and control a large number of 
operations which are required for normal function and the ways in 
which cellular systems fail in disease. 
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Figure 9 | (a): The whole row index of each column of matrix (72* ® 1 ) Hj obtained in system (24). (b): The whole row index of each column of matrix 
(724®lj4)Hj. Each point corresponds to the row index of each column, which implies the position of 1. 
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Figure 10 | (a): Time evolution oftheTBNs (24) with initial state sequences Xi(-l) - (0,0),X2(-1) ^ (1, l),Xi{0) - (1,0),X2(0) - (0, 1). The red solid line 
presents the state evolution of Xi, while the green solid line presents the state evolution ofx2. (b): Total synchronization error E{t) of the TBNs (24). Fig. (a) 
shows states x\ and , states x^ and xl can not reach synchronization till time to 1 00, while Fig. (b) shows the total synchronization error fluctuates between 
value 2 and value 0. Both Fig. (a) and Fig. (b) imply that the coupled TBNs can not be synchronized. 
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